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Abstract 

We present preliminary measurements of branching fractions and charge asymmetries for the B 
meson decays B — ► rf'>K*, B — ► r/^p, and B + — ► r]'ir + . The data were recorded with the BABAR 
detector at PEP-II and correspond to 89 x 10 6 BB pairs produced in e + e~ annihilation through 
the T(4S) resonance. We find the branching fractions B{B° -> rjK*°) = (19.0l|? ± 1.3) x 10~ 6 , 
B(B + -► r/K*+) = (25.7^ + 1.8 )xKT 6 , /3(B+ -> r/p+) = (10.5±| j±1.3) x 10" 6 , B{B + -► r/'p + ) = 
(14.0j£J ± 1.9) x 10~ 6 (< 22 x 10" 6 with 90% confidence), and B(B + -»■ r/7r + ) = (2.8]t^ ± 0.3) x 
10~ 6 (< 4.5 x 10~ 6 ). We also set 90% CL upper limits of B{B° -► r/if*°) < 6.4x 10~ 6 and B(B + -> 
r}'K* + ) < 12 x 10~ 6 . The time-integrated charge asymmetries are A c h{riK* ) = +0.03 + 0.11 + 0.02, 
A ch {r,K* + ) = +0.15 ± 0.14 ± 0.02, and = +0.06 ± 0.29 ± 0.02. 
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1 Introduction 



We report the results of searches for B decays to the charmless final states 1 rf^K* , rj^'p , and 
r]'7r + . For decays that are self tagging with respect to the b or b flavor we also measure the direct 
CP-violating time-integrated charge asymmetry, A c h = (r _ —T + )/(T~ + T + ). For charged B 
decays T^ 1 = T(B ± — > rj^h^), and for B° — > 77.fr* with if* — > if + 7r~ the sign on T matches that 
of the tertiary kaon. 




(c) (d) 

Figure 1: Feynman diagrams for the decays B — ► (rj,r]')(K, K'* ) + . The neutral decays are similar 
except that the spectator quark becomes a d, and the tree diagram is internal. 

Interest in B decays to rj or 77' final states intensified in 1997 with the CLEO observation of the 
decay B — * vj ' K [1]. It had been pointed out by Lipkin six years earlier [2] that interference between 
two penguin diagrams (see Fig. la and lb) and the known 77/7/ mixing angle conspire to greatly 
enhance B — > r/'K and suppress B — > r\K. Because the vector K* has the opposite parity from the 
kaon, the situation is reversed for the B — > 77' if* and -B — > r]K* decays. The general features of 
this picture have already been verified by previous measurements and limits. However the details 
and possible contributions of flavor-singlet diagrams (c./. Fig. Id) can only be tested with the 
measurement of the branching fractions of all four (j],if)(K,K*) decays; the branching fraction of 
the B — > r]'K* decay is expected to be particularly sensitive to a flavor-singlet component [3, 4]. 
In any case the tree diagrams (Fig. lc) are CKM suppressed. The results described in this paper 

1 Except as noted explicitly, we use a particle name to denote either member of a charge conjugate pair. 
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complete the measurement of all four decays with a BABAR dataset of 89 million BB decays [5, 6]. 

The situation for the decays rj^ 1 ' p and t]'tt + is different. These decays are expected to be 
dominated by tree diagrams (see Fig. 2c and 2d) since it is now the penguin diagrams (Fig. 2a 
and 2b) that are suppressed. Since the internal tree diagram (Fig. 2d) is color suppressed, any 
interference effect is diluted for these decays. Branching fractions for these decays are generally 
expected to be in the range (1-10) xl(T 6 [7, 8, 9, 10]. 




»> »» 

u u u u 



(c) (d) 

Figure 2: Feynman diagrams for the decays B + — > rjp + , B + — > r]'7r + , and B + — > rj'p + . 

The charge asymmetry A c h for most of these decays is expected to be small [7, 11]. However, 
for the decay B — ► r/'n and B — > rj'K*, the penguin and tree diagrams are expected to be of similar 
magnitude, which allows potentially large charge asymmetries [4, 9, 10, 12]. 

The current knowledge of these decays comes from published measurements from CLEO [13] 
and BABAR [14], and conference results from BABAR [15] and Belle [16]. Tables 1 and 2 summarize 
these previous results. 

2 Detector and Data 

The results presented in this paper are based on data collected with the BABAR detector [17] at 
the PEP-II asymmetric e + e~ collider [18] located at the Stanford Linear Accelerator Center. An 
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Table 1: Summary of branching fraction results for B decays to rj mesons from CLEO [13], previous 
BABAR [15] measurements, Belle [16], and the present analysis. The results for all fits are given 
as well as a 90% CL upper limit if given in the cited reference. The overall yields and efficiencies 
(e) are given as the sum of yields and efficiencies from the daughter particle decay channels. Rows 
marked (*) refer to measurements that are superseded by the results of the present paper. 
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# BB (10 6 ) 


Fit B(W' 6 ) 


UL B(W' 6 ) 


Signif. (er) 


Signal yield 
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± 1.6 
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15.6 
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BABAR* 
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19.1+ 


6.3 
5.4 
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5.4 


Zi - ' -6.1 
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Belle 


23 


21.2113 


±2.0 
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5.1 


22.1 
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This rpsnlt 

-L lllu A v_' i-> LAX U 
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19.011? 


± 1.3 




15 


157 


9.2 






26.4ll:« 












CLEO 


10 


±3.3 




4.8 


19.2 


7.0 


BABAR* 


23 


22.11ft 1 


< 34 


3.2 


14 ±7 


3.2 


Belle 


23 






< 50 


2.8 


13.8 


2.5 


This result 


89 


nr 7+3.8 
Zd - ' -3.6 


± 1.8 




12 


113 


4.7 


VP + 
















CLEO 


10 


4-811:1 


< 15 


1.3 


5.0 


11.9 


Belle 


23 






< 6.8 


0.0 


0.0 


3.9 


This result 


89 


10.51^ 


± 1.3 




4.8 


65 


7.0 



integrated luminosity of 82 fb _1 , corresponding to 89 million BB pairs, was recorded at the T(45) 
resonance ( "on-resonance" , center-of-mass energy yfs = 10.58 GeV). An additional 9.6 fb _1 were 
taken about 40 MeV below this energy ("off-resonance") for the study of continuum backgrounds 
in which a light or charm quark pair is produced instead of an T(4S). 

The asymmetric beam configuration in the laboratory frame provides a boost of (3j = 0.56 to the 
T(4S). Charged particles (tracks) are detected and their momenta measured by the combination 
of a silicon vertex tracker (SVT), consisting of five layers of double-sided detectors, and a 40- 
layer central drift chamber, both operating in the 1.5-T magnetic field of a solenoid. Photons and 
electrons (neutral clusters) are detected by a CsI(Tl) electromagnetic calorimeter (EMC). 

Charged-particle identification (PID) is provided by the average energy loss (dE/dx) in the 
tracking devices and by an internally reflecting ring-imaging Cherenkov detector (DIRC) covering 
the central region. 

3 Event Selection 

We reconstruct the rj mesons in both of the dominant final states rj — > 77 (t/ 77 ) and 77 — > 7r + 7r~7r° 
(%7r)- For the rf, we also reconstruct two final states: rf — > p°7 (rf^y) and rf — > r/7r + 7r~ (^^), 
with 77 — > 77 (except in the rj^^^^K* mode, where we include also rj — > 7r + 7r~7r°). The K*° 

is reconstructed as K+tt' (K^° +7r _), and K* + as either K+ir° (Kp~ + ) or K°ir + (K*± + ), with 
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Table 2: Summary of branching fraction results for B decays to rf mesons from CLEO [13], previous 
BABAR [15, 14] measurements, Belle [16], and the present analysis. The results for all fits are given 
as well as a 90% CL upper limit if given in the cited reference. The overall yields and efficiencies 
(e) are given as the sum of yields and efficiencies from the daughter particle decay channels. Rows 
marked (*) refer to measurements that are superseded by the results of the present paper. 
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Kg — > 7r + 7r _ . The p + is reconstructed as 7r + 7r°. 

Monte Carlo (MC) simulations [19] of the signal decay modes and of continuum and BB back- 
grounds are used to establish the event selection criteria. The selection is designed to achieve high 
efficiency and retain sidebands sufficient to characterize the background for subsequent fitting. 
Photons must have energy exceeding a threshold dependent on the combinatorial background of 
the specific mode: E 1 > 30 MeV for the two photons used to reconstruct the it in 77 — > 7r + 7r~7r° 
candidates, E 1 > 100 MeV for rj — > 77, and B 7 > 200 MeV for the photon in 77' . Additionally, in 
order to reject background from B — ► X*7 for the r/K* and r\p analyses, we require that the cosine 
of the center of mass decay angle for r/ 77 daughters, relative to the flight direction of the 77, have 
an absolute value of less than 0.86. 

We select resonance candidates with the following requirements on the invariant mass (in MeV) 
of their final states: 910 < < 1000 for rf and rf^^, 520 < < 570 for 773^, 755 < mx* < 
1035, and 470 < m p + < 1070. For the rj'^^^^K* channel, we tighten the K* mass range to 
792 < m K *o < 992. Since these quantities are observables input to a maximum likelihood fit, 
the criteria are loose. Additional states are selected with 2-3 sigma cuts: 120 < m n o < 150, 
490 < m v < 600 for r/ 77 , tightened to 510 < m v < 580 in the ry^ 7r7r p + channel to reduce the amount 
of continuum background in the sample. For Kg — > 7r + 7r _ candidates we require 488 < rriK s < 508. 

In modes with a or a p + , we also require that the cosine H of the helicity angle (the 

vector meson's rest decay angle with respect to its flight direction) be greater than —0.5. For 
/ 3w \ ww K*°, we require H > —0.9. For Kg candidates we require that the lifetime significance 
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(r/<7 T ) be > 3. 

We make several particle identification (PID) requirements to ensure the identity of the signal 
pions and kaons. Tracks in resonance candidates (excluding the Kg) must have DIRC, dE/dx, 
and EMC responses consistent with the expected particle type. For the B + — ► t]'tt + decay, we 
require that the prompt charged track have an associated DIRC Cherenkov angle between — 3.5 a 
and +3.5cr from the expected value for either a pion or a kaon. The Cherenkov angle is used in a 
subsequent fit to distinguish 7]'tt + from rf'K + . 

A B meson candidate is characterized kinematically by the energy-substituted mass mES = 
\J (-^s + p • pb) 2 /Eq — p B and energy difference AE = E* B — \\J~s, where the subscripts and B 
refer to the initial T(4S) and to the B candidate, respectively, and the asterisk denotes the T(4S) 
frame. The mode-dependent resolutions on these quantities measured for signal events average 
about 30 MeV for AE and 3.0 MeV for m ES . We require \AE\ < 0.2 GeV and 5.2 < m ES < 
5.29 GeV. 

3.1 Tau, QED, and continuum background 

To discriminate against tau-pair and two-photon background we require that each event contain 
at least three charged tracks, including at least one from the recoil B meson in addition to those 
required to complete the reconstructed B candidate. 

To reject continuum background, we make use of the angle 9t between the thrust axis of the 
B candidate and that of the rest of the tracks and neutral clusters in the event, calculated in the 
T(4s) frame. The distribution of cosOt is sharply peaked near ±1 for combinations drawn from 
jet-like qq pairs and is nearly uniform for the isotropic B meson decays; we tune the cos #t cut for 
each channel to optimize our sensitivity to signal in the presence of continuum background. The 
resulting requirements are | cos6>t| < 0.8 for r/ 77 p + and , n' J] t 3n \ nn K* , | cos^tI < 0.75 for rj'^K* and 
rj'K* + , | cos 6t\ < 0.65 for 7^ 7 7r + and rj'p + , and | cos 6t\ < 0.9 for all of the other decay sequences. 
The average candidate multiplicity is about 1.1 to 1.2 per event. In events with more than one 
candidate we accept the one with r\ or rf mass closest to the nominal value. 

The remaining continuum background dominates the samples and is modeled from sideband 
data for the maximum likelihood fits described in Section 4. 

3.2 BB background 

We use Monte Carlo simulations of B°B° and B + B~ pair production and decay to look for possible 
BB backgrounds. Most BB backgrounds in these analyses come from a handful of charmless 
decays with similar final states. To study these backgrounds, we perform our preselection cuts 
on high statistics MC samples of the prominent BB backgrounds. Using our resulting selection 
efficiency, along with measured (or predicted) branching fractions for the charmless decays in ques- 
tion, we estimate the number of BB background events expected to enter our samples from each 
exclusive charmless decay. From simulated experiments where we embed the expected number of 
BB background events from the MC, we determine whether these backgrounds are large enough 
to warrant an extra component in the likelihood fit, as described in Section 4.1. 

From these studies we find no evidence for significant BB backgrounds in either the r] — > 
7T + 7T~7r decay chains or the 7/^ 7r7r 7r + decay. For the other channels, however, we see evidence for 
potential BB backgrounds, normally involving other charmless B decays with a high-momentum 
r\ or rf '. Such backgrounds typically are at a level of one or two events and we include a BB 
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component in the likelihood fits to discriminate between these backgrounds and the signal. For 
those channels where we include a BB component in the fit, we generate composite BB background 
samples, typically containing events from about five exclusive charmless decay chains. 



4 Maximum Likelihood Fit 

We use an unbinned, multivariate maximum likelihood fit to extract signal yields for our modes. A 
sub-sample of events to fit for each decay channel is selected as described in Section 3. The sample 
sizes for the decay chains reported here range from 700 to 30000 events. 



4.1 Likelihood Function 

The likelihood function incorporates a number of observables to distinguish signal from the large 
number of background events retained by the sample selection. We describe the B decay kinematics 
with two variables: AE and uies- We also include the mass of the primary resonance (m v or m v i) 
and a Fisher discriminant which describes energy flow in the event. The Fisher discriminant 
combines four variables: the angles with respect to the beam axis of the B momentum and B 
thrust axis (in the Y(4S) frame), and the zeroth and second angular moments Lq^ of the energy 
flow about the B thrust axis. The moments are defined by Lj = J2iPi x |cos#i| J , where 6i is the 
angle with respect to the B thrust axis of track or neutral cluster i, pi is its momentum, and the 
sum excludes the B candidate. 

In addition to these four variables, we include the following observables relevant to particular 
decay chains. For modes with a K* or p, we include in the fit the mass and helicity angle of the 
vector meson with a two-dimensional probability distribution function (PDF) accounting for the 
different shapes of the helicity-angle distributions for true K* or p mesons in the background and 
for combinatoric background. To separate 7/71"+ from 77' ' K + we include in the PDF for the pion 
(kaon) signal the observable S w (Sk), the Cherenkov angle residual with respect to the expected 
angle for pions (kaons) normalized by the measurement error. 

As measured correlations among these observables in the selected data are small, we take the 
PDFs for each event i to be a product of the PDFs for the separate observables. We define hy- 
potheses j, where j can be signal, continuum background, or (where appropriate) BB background. 
The PDF for the rf'^ p and rf'^K* analyses is given by 

V) = V 3 (m ES l ) ■ PjiAE?) ■ V 3 \F) ■ Pjim^) ■ Pj(m* K . /p ,7C). (1) 
For the r]'ir + analysis, j separately indexes pion and kaon components: 

V) = V 3 (mES l ) ■ Vj{U?i) ■ Vj{F) ■ ^K(o) • Vj(Si). (2) 
The likelihood function for each decay mode is 

* 3 

where Yj is the yield of events of hypothesis j to be found by the fitter, and N is the number 
of events in the sample. The first factor takes into account the Poisson fluctuations in the total 
number of events. 
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4.2 Signal and Background Parameterization 

We determine the PDFs for signal from MC distributions in each observable. The PDFs for BB 
background (where appropriate) arise from fitting the composite BB MC sample, discussed in 
Section 3.2. For the continuum background we establish the functional forms and initial parameter 
values of the PDFs with data from sidebands in tties or AE. We then refine the main background 
parameters (excluding resonance mass central values and widths) by allowing them to float in the 
final fit. 

The distributions in the resonance mass(es), ttt-es, an d AE for signal, are parameterized as 
Gaussian functions, with a second Gaussian as required for good fits to these samples. Slowly 
varying distributions (combinatoric background under the resonance mass and AE peaks) are 
parameterized by linear or quadratic functions. We find that two Gaussians also describe the AE 
shape in BB background, with one Gaussian typically peaking low (for feed-down) and the other 
centered high (for feed- up). The combinatoric background in toes is described by a phase-space- 
motivated empirical (ARGUS) function [20]. For BB background, we use an ARGUS function 
plus a Gaussian to parameterize toes- We model the T distribution using a Gaussian function with 
different widths above and below the mean; in background, we include a second Gaussian or a linear 
contribution to account for outlying events. The additional terms ensure that the background is 
well modeled in the low-side tail region where signal is found. 

The PDFs for the helicity angle variable H are given by polynomials (quadratic in signal and 
continuum background, quartic in BB background). This function is multiplied by a Fermi-Dirac 
threshold function where needed to fit a drop in detector efficiency near TL = +1 or — 1 where 
the energy of one of the resonance daughters is low. As there is a true resonance component in 
continuum background, we allow for different shapes in TL for those events corresponding to real 
K* or p resonances as opposed to the combinatoric continuum. 

The PDFs T-^SV,^) are determined with a sample of D*-tagged D° decays, and parameterized 
as double Gaussians. Consistent with the calibrations, the kaon PDFs for S and AE are copies of 
those for pions evaluated at the corresponding displaced value of the observable 2 . 

We check the simulation on which we rely for signal PDFs by comparing with large data control 
samples. For to.es and AE we use the decays B~ — > ir~D° and B~ — > p~D° with D° — > K~7r + 7r°, 
which have similar topology to the modes under study. For the resonance masses we use inclusive 
resonance production in data. We adjust the means and widths of PDF parameterizations based 
on these control samples. 

5 Fit Results 

By generating (from PDF shapes) and fitting simulated samples of signal and background, we verify 
that our fitting procedure is functioning properly. We find that the minimum — In C value in the 
on-resonance sample lies well within the — ln£ distribution from these simulated samples. 

The efficiency is obtained from the fraction of signal MC events passing the selection, adjusted 
for any bias in the likelihood fit. This bias is determined from fits to simulated samples, each equal 
in size to the data and containing a known number of signal (and BB, where appropriate) MC 
events combined with events generated from the continuum background PDFs. The biases we find 

2 To facilitate simulations and PDF projections we actually work with the transformed pairs (S^jSS = S n — Sk) 
and (AE„,8AE = AE K - AE*) 
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depend on the mode, but are a few percent where the yield is substantial, or at most 2-3 events in 
the limit of vanishing yield. 

Table 3: Yields, efficiencies, branching fractions, and charge asymmetries for B decays to states with 
an rj meson, measured with our sample of 89 million BB pairs. The overall yields and efficiencies 
(e) of the combined mode are given as the sum of yields and efficiencies from the various daughter 
particle decay channels, and for these we incorporate the systematic errors. 



Mode 


Fit £(10~ 6 ) Signif. (a) Signal yield e (%) A ch 


B — > r]K 


19.U_2.! ± 1.3 15 157 9.2 +0.03 ± 0.11 ± 0.02 
21 ±3 13.1 125 ±16 6.6 ±0.12 ±0.12 
14 ±4 6.8 32 ±9 2.6 -0.39 ± 0.25 


b+ n K* + 


25.7±f;§ ±1.8 12 113 4.7 ±0.15 ± 0.14 ± 0.02 
23 ±6 5.7 46 ±12 2.2 ±0.04 ± 0.24 

32 ±9 6.2 27 ±8 0.9 +0.46±g;|| 

27 ±8 6.0 30 ±9 1.1 — 0.06±g;|f 
20 ±10 4.5 10 ± 5 0.5 +0.39±g;|J 


B+ np + 

V3ttP + 


10.5l|g± 1.3 4.8 65 7.0 ±0.06 ± 0.29 ± 0.02 

10 ±4 3.8 42 ±17 4.7 ±0.46±^ 

11 ±6 3.0 23 ± 12 2.3 -0.56^;^ 



In Tables 3 and 4 we show the results of the fits to the on-resonance data. Shown for each 
decay mode are the measured branching fraction and the 90% confidence level upper limit (where 
appropriate). We report the statistical significance for the individual decay chains and display 
the significance including systematics for the combined result in each channel. The statistical 
significance is taken as the square root of the difference between the value of — 21n£ for zero signal 
and the value at its minimum. The tables also include the fit signal yield, where the statistical 
error on the number of events is taken as the change in the central value when the quantity 
— 21n£ changes by one unit. The 90% C.L. upper limit is taken as the solution B to the condition 
Jo C(b)db/ Jo 00 C{b)db = 0.9. The overall efficiency is given by the product of the Monte Carlo 
efficiency, efficiency corrections, and branching fraction products for the daughter particle decay 
sequences. The number of produced B mesons is computed with the assumption of equal production 
rates of charged and neutral pairs. The final column in Table 3 gives the charge asymmetry (A c h)- 

In Figures 3 and 4, we show projections of vies and Ai? for the nK* and i]p + modes, respectively. 
We make these plots by selecting events with signal likelihood (computed without the variable shown 
in the figure) exceeding a mode-dependent threshold that optimizes the expected sensitivity. The 
selection retains a fraction of the signal yield averaging about 70% across the decay sequences. 

6 Systematic Uncertainties 

Floating background parameters in the fit allows us to incorporate into the overall statistical error 
most of the systematic errors on yields that arise from uncertainties in the values of the PDF 
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Figure 3: Projections of the B candidate ttt-es and AE for r]K*° (a, b) and r]K* + (c, d). Points 
with errors represent data, solid curves the full fit functions, and dashed curves the background 
functions. These plots are made with a cut on the signal likelihood and thus do not show all events 
in the data samples. 
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Table 4: Yields, efficiencies, and branching fractions for B decays to states with an rf meson, 
measured with our sample of 89 million BB pairs. We quote 90% CL upper limits if the significance 
of the measured yield is less than 4 sigma. The overall yields and efficiencies (e) of the combined 
mode are given as the sum of yields and efficiencies from the various daughter particle decay 
channels, and for these we incorporate the systematic errors. 



Mode 


Fit £(KT 6 ) UL B(10- e ) Signif. (a) Signal yield e (%) 


H f 


2.8±i g ± 0.3 < 4.5 3.4 13 9.9 
4±2 3.9 17 ±7 4.9 
-1±3 -4 ±10 5.0 


B° tj'K* 
rl K*° 

'?7(77)7r 7T 

rj'/o n K*° 


3.2+i;| ± 0.9 < 6.4 2.2 22 6.6 
-2 ±2 -4 ±4 2.3 
13 ±6 2.8 11 ±5 1.0 
5±3 1.9 15 ±10 3.4 


B+ -► ri'K* + 
rl K* + 

rl K* + 

r lfn iy K s n+ 
rl K* + 
rl K* + 


6.1+|;| ± 1.2 < 12 2.4 15.4 3.0 
-12 ±5 -8 ±4 0.8 

16 ±9 2.6 16 ±9 1.1 

7±7 1.8 3±3 0.5 
8 ± 13 0.7 5±7 0.6 


b+ n' P + 

Vr/mrP 
V'fryP + 


u.otli ± 19 < 22 3 - 8 69 41 
11 ±6 3.1 17 ±8 1.8 
25 ± 12 2.8 52 ± 23 2.3 



parameters. We determine the sensitivity to parameters of the signal PDF components by varying 
these within their uncertainties in all channels except for r}'ir + . In the r/7r + decay, we float the 
primary signal parameters and determine the effect on signal yield; this is possible because of the 
simultaneous fit with f]'K + , which has a large branching fraction. This is the only systematic error 
on the fit yield; the other systematics apply to either the efficiency or the number of BB events. 

The uncertainty in our knowledge of the efficiency is found to be 0.8N t %, 2.5iV 7 %, and 3% for a 
Kg decay, where N t and iV 7 are the number of signal tracks and photons, respectively. We estimate 
the uncertainty in the number of produced BB pairs to be 1.1%. The estimated error on systematic 
bias from the fitter itself (1-12%) comes from fits of simulated samples with varying signal and 
background populations. Published world averages [21] provide the B daughter branching fraction 
uncertainties. We account for systematic effects in cos#t (1 — 3.5%, depending upon how tight a 
cut is made), in the PID requirements (2%), and from MC statistics (1.2%). For decay channels 
where a BB component is included in the fit, we assign an additional systematic determined by 
varying the BB yield within its errors. We quote the BB background systematic to be one half 
the measured change in the signal yield when the number of fit BB events is changed by one 
standard deviation. For the ri'p + decay, where the fit BB background yield is about five times 
larger than other decays, we also vary the BB background model to account for different amounts 
of the expected individual backgrounds. 
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Figure 4: Projections of the B candidate toes ( a ) an d AE (b) for B + — > Points with errors 
represent data, solid curves the full fit functions, and dashed curves the background functions. 
These plots are made with a cut on the signal likelihood and thus do not show all events in the 
data samples. 



A study of the charge asymmetry as a function of momentum for all tracks in hadronic events 
bounds the tracking efficiency component of charge-asymmetry bias to be less than 1%. _D*-tagged 
D — > Kit and B samples provide additional crosschecks that the bias is small. We assign a 
systematic uncertainty for A c h of 2% based on the tracking study and a small PID contribution 
determined from the D* studies. 

We keep track of which systematic error contributions are (un) correlated between the several 
measurements with different secondary decay modes of the same primary decay for use in obtaining 
their average as the final branching fraction result. 



7 Combined Results 

To obtain the final results we combine the branching fraction and charge asymmetry measurements 
from the individual daughter decay chains. The joint likelihood is given by the product, or equiva- 
lently — 21n£ is given by the sum, of contributions from the submodes. The statistical contribution 
comes directly from the likelihood fit, which reflects the non-Gaussian uncertainty associated with 
small statistics. Before combining we convolve each statistical — 21n£ with a Gaussian function 
representing the part of the systematic error that is uncorrelated among the submodes. We show 
the resulting distributions in Fig. 5 for our measurements of previously unseen decays that have 
significance greater than 3 sigma. The corresponding distributions without systematics give the 
combined statistical errors, and these in conjunction with the solid curves in Fig. 5 and the corre- 
lated systematics give the total systematic errors. 

The resulting branching fractions and charge asymmetries are included in Tables 3 and 4, where 
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Figure 5: Distributions of —2 In £ convolved with uncorrelated systematic errors for branching frac- 
tion measurements. Solid curves represent the result of combining channels. The frames present (a) 
ijp + , with r? 77/ o + (dashed) and r]3 n p Jr (dotted); (b) r/p + , with rj' vvn p + (dashed) and rj' n p + (dotted); 
and (c) i]'tt + , with rj n ir + (dashed) and ir + (dotted). 
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the significance given includes systematics. 



8 Conclusion 

We report preliminary measurements of branching fractions and A c h for B meson decays to r/ 
or rj with a K*, p + , or 7r + . We find signals with statistical significance exceeding four standard 
deviations in all r] channels. The decay B + — > r/p + has not been seen previously. We also have 
evidence for B + — ► r/'p + (with significance 3.8<t) and B + — > 7/71-+ (with significance 3.4cr). The 
observed values in the rj channels are 

B(B° -► r/K* ) = (19.0±£? ± 1.3) x 10~ 6 , 
B{B + ^r]K* + ) = (25.7±1;| ± 1.8 ) x HT 6 , 
B(B + -► = (10.51^ ± 1.3) x 10~ 6 . 

For the rj channels we find 

= (2.8±i;g ±0.3) x 10~ 6 (< 4.5 x 10~ 6 ), 

< 6.4 x 10~ 6 , 

< 12 x 10~ 6 , 

= (14.0±2:e ± L9 ) x 10 ~ 6 (< 22 x 10 ~ 6 )> 

where the upper limits are taken at 90% CL. These results supersede the previous BABAR measure- 
ments [14, 15]. They represent substantial improvements over all previous measurements, as can 
be seen from Tables 1 and 2. The branching fraction limit for B + — > t]'tt + is nearly three times 
more restrictive than previous measurements. The measurement for B — > r/'K* is not yet precise 
enough to determine whether a flavor singlet component is present for this decay, though we do 
restrict the size of this contribution. 

For the modes with significant signals, we measure the charge asymmetries 

A ch (r]K*°) = +0.03 ±0.11 ±0.02, 
A ch (r]K* + ) = +0.15 ±0.14 ±0.02, 
A ch {rip + ) = +0.06 ±0.29 ±0.02. 

These charge asymmetry results are in agreement with the theoretical expectations discussed in 
Section 1 and rule out substantial portions of the physical region. 
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